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Optical Weights and Waterfalls in Doped Charge Transfer Insulators: an 

LDA+DMFTStudy of LSCO 

Cedric Weber/ Kristjan Haule/ and Gabriel Kotliar-'^ 
^Department of Physics, Rutgers University, Piscataway, NJ 08854, USA 

We use the Local Density Approximation in combination with the Dynamical Mean Field Theory 
to investigate intermediate energy properties of the copper oxides. We identify coherent and inco- 
herent spectral features that results from doping a charge transfer insulator, namely quasiparticles, 
Zhang-Rice singlet band, and the upper and lower Hubbard bands. Angle resolving these features, 
we identify a waterfall like feature, between the quasiparticle part and the incoherent part of the 
Zhang-Rice band. We investigate the assymetry between particle and hole doping. On the hole 
doped side, there is a very rapid transfer of spectral weight upon doping in the one particle spec- 
tra. The optical spectral weight increases superlinearly on the hole doped side in agreement with 
experiments. 



Since their discovery, the high temperature supercon- 
ductors continue to be a subject of intensive investiga- 
tions. It is widely believed that strong correlations and 
many body effects are responsible for many of the pecu- 
liar properties of these materials. However, after many 
years of intensive studies, a comprehensive understanding 
of their electronic structure, even at intermediate energy 
scales, is lacking. The question of wether the strengh 
of the on site correlations is sufficient to open the gap 
at half filling (Mott picture) or alternatively magnetism 
is required (Slater picture) is actively debated [T]. A 
vertical feature in the photoemission spectral intensity, 
commonly known as a waterfall, has recently been found 
in many compounds [2", 3 , and it has been assigned to 
several conflicting microsopic origins [H [51 [Hj • 
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It is generally accepted that these materials contain 
copper-oxygen layers as a common element. Stoichiomet- 
ric compounds such as La2Cu04 , are antiferromagnet- 
ically ordered charge transfer insulators in the Zaanen- 
Sawatsky Allen [3 [8] classification scheme. Hence a min- 
imal model of these materials contains two oxygen and 
one copper orbital per unit cell [9j Ji)^ . When this model 
is doped with electrons, the low energy physics is be- 
lieved to be closer to that of the Hubbard model. On the 
hole doped side, this reduction constructed by Zhang and 
Rice [iTl [12] can be achieved in the limit that the charge 
transfer energy is much larger than the copper oxygen 
hybridization. However, there is no agreement on the 
precise region ( in energy, and parameter space ) where 
the reduction to a one band model is accurate. 

In this letter, we use a realistic theoretical approach, 
i.e., the Local Density Approximation combined with 
the Dynamical Mean Field Theory (LDA-hDMFT) [13], 
to study a typical cuprate compound La2_a;Sra;Cu04 . 
The LDA calculation was done by PWSCF package [Tij . 
which employs a plane-wave basis set and ultrasoft pseu- 
dopotentials [T5]. Downfolding to a three band model, 
containing copper dx2-y2 and two oxygen p^- orbitals was 
performed by the maximally localized Wannier functions 
(MLWF) method HH [T7]. We solve this model using 



FIG. 1: (Colors online) Doping S (left scale) versus the chem- 
ical potential fi. The parent compound {5 = 0) is a charge 
transfer insulator with a gap Apd ~ 1.8eV (green region). 
The number of electrons in the di,2-i/2 orbital (n^) and in the 
Per orbitals (up) are also shown (right scale). 



Dynamical Mean Field Theory with an exact impurity 
solver [18]. 

We address the following questions. What are the spec- 
tral features that result from doping a charge transfer 
insulator? How are these spectral features distributed in 
momentum space? How does the quasiparticle residue 
and the optical weight, induced by doping the charge 
transfer insulator, vary with doping? We address the 
similarities and differences between our results and those 
that follow from a Hubbard model description. 

Downfolding the LDA band structure of La2Cu04 re- 
sults in the following three band Hamiltonian: 



ijcr,(a,fl)e{pj:,py.d^2-y2) 
iaae(p^,Py) 
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where i and j label the Cu02 unit cells of the lattice, and 



,a/3 



tj are the hopping matrix elements, e^ and tp are the 
on-site energies of the d and p orbitals, respectively. 
The charge transfer energy, Ap^j, between the copper 



and oxygen plays the role of an effective onsite repulsion 
t/ in a Hubbard model picture, as seen for example in 
slave bosons mean-field studies [19]. 

The LDA downfolding procedure results in Cd — ^p — 
2.7^eV . To this hamiltonian, we add the onsite Coulomb 
repulsion U on the dx2~y2 orbital 



"Hu — Ud 2_^ nid^fi-idi 



(2) 



where the value of Ud = SeV . The LDA+DMFT method, 
accounts for the correlations which are included in both 
LDA and DMFT by a double counting correction to the 
d-orbital energy, which we take to be fixed at Edc — 
3.12eT/ for all dopings. We neglect the oxygen Coulomb 
repulsion Up. Within single site DMFT, the copper- 
oxygen repulsion Vdp is treated at the Hartree level, and 
hence it is included in the orbital energies of the p and 
the d orbitals. 

The Green function of the three band model is given 
by: 



Gk(ia;„) = {iun + /i - Hk - S(itJ„))" 



(3) 



where H^ is the Fourier transform of the Jit in Eq. M 
and is a 3 X 3 matrix. S is the self-energy matrix being 
nonzero only in the d orbital. 

The self energy in Eq. Q is obtained by solving an 
Anderson impurity model subject to the DMFT self- 
consistency condition: 



iu) — Eimp — S(icL;) — A{iuj) 



1 



E Gi'i^u;), (4) 
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where the sum runs on the first Brillouin Zone (BZ). In 
this work we use the continuous time quantum Monte 
Carlo impurity solver algorithm [T8j [21] , which gives the 
self-energy functional S[i?imp, A] within the residual sta- 
tistical error bars of the Monte Carlo algorithm [32j Real 
frequency resolved quantities were obtained by analytic 
continuation (more details will be given elsewhere [H]) of 
the observables on the imaginary axis. We have crossed 
checked the analytic continuation using the OCA real 
frequency solver [23j . 

Equations ^ and Q were studied previously in 
Refs. [m US). When Ud is large, there is a metal to 
charge transfer insulator transition at integer filling, as a 
function of the charge transfer energy e^ — Cp. For the set 
of realistic parameters considered in this letter, we find 
that the parent compound is a charge transfer insulator 
(see Fig. [T|, with a hole density per Cu atom « 85% 
and a charge gap close to l.SeV. This places this mate- 
rial slightly above (but not very far) from the metal to 
charge transfer insulator transition point. This is demon- 
strated in Fig. Ill where the number of particles (i5) as a 
function of chemical potential (/i) exhibits a plateau in 
the interval fi = 14:.2eV—\&eV . This figure also display 
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FIG. 2: (Color online) a) Left panel: Spectral function 
A{]<i,uS) for the hole doped compound {5 — 0.24). The black 
lines correspond to the LDA bands. Right panel: The local 
spectral functions of the dx2-y2 band (dashed line) and of the 
Per band (full line). We marked the positions of the lower Hub- 
bard band (L.H.B.), coherent part of the Zhang-Rice singlet 
called the quasi-particle peak (QU.P), the incoherent part of 
the Zhang-Rice singlet band (Z.R.S.) and the upper Hubbard 
band (U.H.B.). b) The low energy part of ^(k, u) from panel 
(a). The dashed line shows the slope of the quasiparticle band, 
which gives the Fermi velocity. The derivative was obtained 
by a quadratic expansion of 'S{uj) at the Fermi level. The 
Fermi velocities in LDA-I-DMFT is vdmft = 2.15(5) AeV in 
very favorable agreement with experiment {vexp = 2.0 AeV 
at doping S = 0.24 [20]). The ratio between the DMFT and 
the band velocity is thus vdmft/vlda = 0.49. c) Fermi sur- 
face and the path along which A(k, uj) is shown in (a) and 
(b) d) The density of states in the electron doped compound. 
A very sharp quasi-particle peak develops close to the upper 
Hubbard band in the electron doped case. 



the partial occupancies, describing the relative distribu- 
tion of the occupancies among copper and oxygen. 

We have computed the spectral functions for both the 
parent and a doped compound. In the parent compound 
(not shown), we found four spectral peaks: i) the upper 
Hubbard band around -1-2.5 eV, ii) a singlet bound state 
of oxygen and copper known as the Zhang-Rice singlet 
just below the Fermi level, iii) the oxygen set of bands 
centered around -7eV and iv) the lower Hubbard band at 
-8eV. By hole doping (see Fig. l2]i-right) , the Zhang-Rice 
band is split into two peaks, the quasi-particle peak and 
the incoherent part which is associated with the lower 
Hubbard band in a single-band Hubbard model. 

Figures [2k-left, and I2I3 show the momentum resolu- 
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FIG. 3: (Color online) a) Optical spectral weight A^e// ob- 
tained by Eq. IT] using various cutoffs A = 1.2, 1.5, l.ieV for 
both the hole doping [5 > 0) and electron doping (5 < 0). The 
experimental data (red circles) are reproduced from Ref. [5B] . 
Long dashed line corresponds to the curve N^jj (5) = 5. Note 
the assymetry between the electron [5 < 0) and hole doped 
{5 > 0) side, b) The quasi-particle weight Z versus doping. 
Inset: double occupancy of the dx2-y2 orbital versus dop- 
ing. For the electron doped compound, the double occupancy 
increases linearly with doping, whereas the hole doping has 
only a weak effect on double occupancy of the dx2-y2 orbital. 
Continuous lines are guide to the eyes. 




tion of these features. The quasiparticle band crosses the 
Fermi level with a Fermi velocity along the nodal direc- 
tion of the order of vpermi = 2.2AeV in good agreements 
with experimental results [20] . Figure jSK also shows a 
spectral gap around -0.3 eV, separating the coherent and 
incoherent parts of the Zhang Rice band, with weak spec- 
tral weight connecting them. This is reminiscent of the 
waterfall feature observed in numerous cuprates. More- 
over, our calculations lead to a very good agreement with 
experiments for both the range of energy where the high 
energy peak is present (~ —leV) and the position in 
the Brillouin zone where the abrubt change is observerd 
(close to (7r/2,7r/2)). While there are many interpreta- 
tions of this feature [3 |3l HI |5l [6] , we find that, a proper 
description of the waterfalls requires a multi-band model 
containing both oxygen and copper orbitals, as pointed 
out in Ref. |4| . The waterfall comes from the spectral gap 
between the incoherent part and the coherent (quasipar- 
ticle) part of the Zhang-Rice band. The latter is not 
present in the single site DMFT description of the par- 
ent compound, and our calculation in paramagnetic state 
of La2Cu04 does not show splitting of the Zhang- Rice 
band into a low energy and a high energy part. We have 
found however, that when the anti- ferromagnetic broken 
symmetry state is considered for a parent compound, the 
splitting of the Zhang- Rice singlet is present, and the low 
energy part of the Zhang-Rice band is touching the Fermi 
level at (7r/2, 7r/2) momentum. 

One of the most notable aspects of cuprate physics, is 
the rapid growth of the optical conductivity inside charge 
gap upon hole doping [26 . This was noted early on by 
Sawatzky and collaborators [22,. In order to tackle this 
issue, we computed the optical conductivity, given by: 



crk 



dx 



f{x - w) - f{x) 



Tr(/5kcr(a; -cj)vk/5kcr(a;)vk) (5) 



Where c is the interlayer distance, the density matrix p 
is defined by 



Pk<T(a;) = — \G\^„{x) - Gk^(x) 



(6) 



FIG. 4: (Color online) Optical conductivity (j{uj) for the 
parent compound [5 = 0) and the hole doped compound 
[5 = 0.24). For the parent compund, we show optics in the 
Neel state. The optical gap in the parent compound is around 
1.8 eV in both the paramagnetic and antiferromagnetic state. 
The structure around ui = 8eV corresponds to transitions be- 
tween the lower and the upper Hubbard band. Upon hole 
doping, additional peak appears close to uj ~ leV for a large 
range of doping 5 = 5 — 30%. It corresponds to transitions 
between the quasi-particle part and incoherent part of the 
Zhang-Rice singlet. 



And the bare vertex is given by Vk — diik/dkx. Our 
results are dispayed in Fig. |4J The optical gap in the 
parent compound is found to be l.SeV. To quantify the 
rate of the redistribution of optical spectral weight, we 
computed the effective electron number per Cu atom de- 
fined by 



</ 



2meV 



a'{uj)d(jj, 



(7) 



where rUe is the free electorn mass, and V is the cell Vol- 
ume containing one formula unit. N^ff is proportional 



to the number of electrons involved in the optical ex- 
citations up to the cutoff A. Our results for -/Ve// are 
displayed in Fig. |3] and compared to experimental data 
taken from Ref. [2B]. 

Notice that the spectral weight, contained in the region 
below the charge transfer gap of the parent compound, 
grows in a superlinear fashion. This is in agreement with 
experiments of Ref. |26j but is incompatible within a rigid 
band picture of doping either a band or a Mott insula- 
tor far from the Mott transition (as for example in the 
Hubbard model well above C/c2)- Thus the optical spec- 
tral weight in the charge transfer insulator grows faster 
with doping than in the Mott insulator, where the linear 
growth of the spectral weight is expected. 

The realistic three-band model within single site 
DMFT approach leads to a considerable asymmetry in 
the electron and hole doped side of the phase diagram. 
The superlinear behavior found in our calculation of the 
three band model is sufficient to explain the experimental 
findings, even though the correlations are strong enough 
to open a charge transfer gap in the parent compound 
even in the paramagnetic state. 

We verified that this interesting particle hole as- 
symetry is present also in the quasiparticle residue Z = 



1- 






and is shown in Fig. 



3 3. This is a one 



electron quantity that can be obtained~directly from the 
imaginary axis data without the need to invoke analytic 
continaution, but still follows the same trend as -/Vg//. 

It is interesting to compare the asymmetry found here 
with the recent results of Ref. "55] on the Anderson lattice 
model. The authors found a much larger particle hole 
assymmetry compared to our results. Most likely this is 
due to a momentum independent hybridization function, 
which makes the Kondo scale in the hole doped side very 
small [55], and hard to reach with conventional QMC. 

Phenomena in the underdopcd region, such as the nor- 
mal state pseudogap and the formation of fermi arcs, are 
not described by single site DMFT and require cluster 
extensions of the formalism [30. However, single site 
DMFT is still expected to describe well intermediate en- 
ergy phenomena, such as the distribution of the spec- 
tral weight upon doping the Mott /charge transfer insu- 
lator, which is the focus of our work. Momentum resolved 
quantities become more accurate in the overdoped region, 
where the self energy is more local. 

For this reason, we have not considered here the de- 
tails of the momentum dependent photoemission spectra 
at low doping, or the details of the frequency dependent 
optical conductivity, which are considereably infiuenced 
by cluster effects. We have rather focused on quantities 
such as Neff, which are less sensitive to cluster correc- 
tions and we rather presented the detailed photoemission 
spectra in the overdoped region, where actually the local 
approximation is valid |30j . 

In conclusion, we have carried out a realistic 



LDA-t-DMFT calculations accounting for the dx2-y2 and 
Per orbitals of the three band model. The electron doped 
side is behaving mainly like the single-band Hubbard 
model. On the other hand, the hole doped side of the 
phase diagram is very different, a) We find a superlinear 
optical spectral weight transfer in agreement with exper- 
iments, b) On the spectral function, we find a spectral 
gap between the quasiparticle band and the incoherent 
Zhang Rice singlet band. This leads to a very abrupt 
and nearly vertical change in the spectral intensity in the 
nodal direction of the spectrum (waterfall) , that ends up 
below 1 eV. 

Our LDA-I-DMFT study support the view that the hole 
doped cuprates are above, but not very far from the metal 
charge transfer insulator transition in agreement with the 
conclusion of earlier slave boson studies [TS] . 

It would be interesting to apply the same methodol- 
ogy to the electron doped cuprates which have a smaller 
charge transfer gap. These materials might have a charge 
transfer energy which lies below the critical value needed 
to sustain a paramagnetic insulator at integer filling, and 
hence might undergo a different metallization process 
upon doping, as suggested in Ref. jJTj. 
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